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Abstract — In this paper we introduce a class of Markov models, termed 
best-effort networks, designed to capture performance indices such as mean 
transfer times in data networks with best-effort service. We introduce the 
so-called min bandwidth sharing policy as a conservative approximation 
to the classical max-min policy. We establish necessary and sufficient er- 
godicity conditions for best-effort networks under the min policy. We then 
resort to the mean field technique of statistical physics to analyze network 
performance deriving fixed point equations for the stationary distribution 
of large symmetrical best-effort networks. A specific instance of such net- 
works is the star-shaped network which constitutes a plausible model of 
a network with an overprovisioned backbone. Numerical and analytical 
study of the equations allows us to state a number of qualitative conclusions 
on the impact of traffic parameters (link loads) and topology parameters 
(route lengths) on mean document transfer time. 

Keywords — best-effort service, max-min fairness, min policy, mean field, 
star-shaped network. 



I. Introduction 

Consider a network handling data flows from several users, 
and assume no quality of service commitments (such as min- 
imum bandwidth allocations) have been made by the network 
to the users. Such a situation has been prevalent in the Internet 
until now, and is likely to remain so for another few years. 

The preferred service model in that situation, known as best 
effort service, consists in allocating a fair proportion of band- 
width to contending users; see, e.g., Bertsekas and Gallager (T). 
There are actually several possible notions of fairness available 
for this bandwidth allocation problem (see, e.g., Mo and Wal- 
rand [2| for a parametric family of fairness criteria covering all 
other notions proposed so far), although the classical notion pro- 
posed in [ 1 1 is the so-called max-min fairness. 

Recent work has led to a relatively good understanding of 
how bandwidth is shared between network users when a given 
congestion control algorithm is used; see, e.g., Massoulie and 
Roberts [3| and references therein. The question of what type 
of fairness is achieved in the current Internet, where Jacobson's 
congestion avoidance algorithm — as implemented in TCP — is 
responsible for congestion control, has been studied in depth by 
Hurley et al. j4j. These studies all assume the number of flows 
remains fixed. 

In comparison, there is little work accounting for the random 
nature of traffic and its impact on user perceived quality of ser- 
vice. Consider for instance the transfer of digital documents 
(Web pages, files, emails,. . . ) using a transport protocol like 
TCP. This constitutes the bulk of Internet traffic today. The 
performance criterion relevant to such transfers is the overall 
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document transfer time. This time is clearly highly dependent 
on the number of ongoing transfers on the links shared by the 
considered connection. This number varies as a random process 
as new connections are established and existing ones terminate 
in a way which depends on how bandwidth is allocated, as well 
as on the underlying traffic parameters. 

In the case of a single bottleneck resource shared perfectly 
fairly, simple traffic assumptions of Poisson arrivals and iden- 
tically and independently distributed document size lead to a 
processor sharing queueing model 0. This fluid flow model 
provides useful results on expected response times as a function 
of the load of an access link or a Web server, for instance. It 
also shows how a form of congestion collapse can occur when 
demand (arrival rate x mean document size) exceeds capacity. 
The processor sharing queue is then no longer ergodic leading 
to unbounded response times. To understand the impact of mul- 
tiple bottlenecks and to investigate the effect of different sharing 
strategies, one would like to dispose of similar analytical results 
for multiple resource systems. 

To the best of our knowledge, the only analytical results avail- 
able so far are in Massoulie and Roberts [5], where the so-called 
linear network topology is investigated. Simulation results for 
the linear network can be found both in [5 1 and in de Veciana et 
al. 0. The main motivation for the present paper is to study the 
performance of best-effort networks with alternative topologies, 
the ultimate objective being the derivation of heuristics enabling 
the performance evaluation of bandwidth sharing in a general 
network. 

In the present paper we report the results of our preliminary 
investigations. These include an analysis of the stability con- 
ditions under which the expected response time remains finite 
in a general network. We also apply mean field techniques to 
evaluate the performance of large symmetrical networks. Nu- 
merical results derived from the model illustrate how response 
times depend on the number of bottleneck links and their uti- 
lization. These results are of some practical interest and aide 
our understanding of the behavior of best effort networks. A 
further significant contribution is the insight provided into the 
inherent difficulty of deriving performance estimates when more 
than one bottleneck limits throughput. 

Section [II] introduces a general class of Markov models for 
best-effort networks which is intended to capture the impact 
of network topology, traffic parameters and bandwidth sharing 
(fairness) criteria on document transfer times. A brief account 
of the results obtained in [5 | is given, and the so-called "min" 
bandwidth allocation is introduced as a conservative approxima- 
tion to max-min fairness. Section iHll then establishes the nec- 
essary and sufficient ergodicity criteria for best-effort networks 
under the "min" policy. Section|lV]introduces the so-called "star 
topology". Its relevance as a model of real networks is discussed 
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Fig. 1 

The linear network 



and a mean field heuristic is proposed. This heuristic is expected 
to be accurate in the asymptotic regime where the number of star 
branches is large. The derived fixed point equations are inves- 
tigated numerically in Section [V] Simulation is used to verify 
the accuracy of the heuristics. Extensions to the star-shaped net- 
work are also considered in Sections [TV] and [V] which notably 
allow an evaluation of the impact of the number of bottlenecks 
on the mean transfer time. 

II. Best-effort networks 

Consider the following network model: a set L of links is 
given, where each link IgI has an associated capacity or band- 
width Q > 0. A set %^ of routes is given, each route being iden- 
tified with a subset of links. Fig.Q]iilustrates the so-called linear 
network: it consists of L links with equal capacity, route which 
crosses each link, and routes r = l,...,L which cross a single 
link. 

To each route r are associated two parameters: X r is the arrival 
rate of new transfer requests along route r, and a, is the aver- 
age document size. We make the following standard simplifying 
assumptions: requests for document transfers along route r ar- 
rive at the instants of a Poisson process with intensity X r , while 
the corresponding document sizes are mutually independent, in- 
dependent of the arrival times, and drawn from an exponential 
distribution of mean a,-. 

These traffic assumptions make the process specifying the 
number of transfers in progress on different routes Markovian 
(see below) and thus greatly simplify analysis. The Poisson ar- 
rivals assumption is not unreasonable in a large network. In 
view of the insensitivity of performance results for an isolated 
link to the exact document size distribution, we do not expect 
divergence of the real distribution from the exponential size as- 
sumption to invalidate the derived conclusions. However, the 
main reason for assuming an exponential distribution is clearly 
one of analytical tractability. 

The network state is summarized by the variables X = {jc r , r £ 
where x r denotes the number of transfers in progress along 
route r. It remains to specify at what speed documents are trans- 
mitted in any given state X in order to turn X into a Markov pro- 
cess with well defined dynamics. Indeed, given the rate ^{X) 
at which documents along route r are transferred when the net- 
work state is X, X is a Markov process with non-zero transition 
rates given by 

J x r — > x r + 1 : rate X r , 

\ x r —>x r — l: rate x r t >r {X)/o r . 

A natural assumption would be to consider that each document 



receives its fair share of bandwidth. For instance, if as in JT] 
fairness is understood as max-min fairness, each transfer along 
route r receives a bandwidth share ^™ m , where 

£^r<Q,V^eA (l) 

r3i 

and for every route r, there is at least one link I G r such that 

£ = Q, and Cr - maxCr- (2) 

These two conditions uniquely determine the bandwidth shares 
^5" m . Having specified the Markov process X, one can then at- 
tempt to study its steady state properties, identifying the condi- 
tions on the load parameters 

Pe = -prJ\^ rCr 

under which it is ergodic and, when it is, determining the sta- 
tionary distribution. Mean transfer times T r along each route r 
can then be computed using Little's law: T r = lEx r /X r . 

It turns out that explicit formulas for steady state distributions 
are typically beyond reach. A notable exception is the linear net- 
work, with bandwidth shares being allocated to realize propor- 
tional rather than max-min fairness; see Q. In order to obtain 
formulas in other cases, one therefore has to resort to asymp- 
totics on various parameters. For instance, for the linear network 
with max-min fair rate sharing, the regime where the arrival rate 
Xq along route goes to zero (essentially, a form of light traffic 
analysis) is considered in [5]; this leads to approximate formu- 
las for 7b. It can be shown, in particular, that 7b increases as 
the logarithm of the number of links L when L increases. This 
is in contrast to the case of proportionally fair sharing where it 
increases linearly in L. 

The main purpose of this paper is to investigate an alterna- 
tive asymptotic regime where it is the network topology which 
evolves. The precise description of this limiting regime will be 
given in Section llVl 

In the following sections, we consider bandwidth allocations 
according to the following "min" policy: given the network state 
X, each transfer along route r receives a bandwidth share i^ 1 ™" 
given by 

r = mm TT> ( 3 ) 

ler X( 

where we have introduced the notation Xf — Y. r J 3l x r l to repre- 
sent the total number of transfers making use of link £. 

It is easy to check that this policy satisfies the capacity con- 
straints (Q]). Moreover, it is sub-optimal with respect to the max- 
min fairness policy, as shown in the next theorem. 

Theorem 1: Under the same initial conditions, the vector 
X mm (t) for the system under the i^ mm allocation policy is 
stochastically smaller than X mm (t), corresponding to the i^ min 
allocation. 

Proof: Assume that, for some f, X mm (t) < X mia (t). Then, 
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with the notation of (0, 



= maxCT 

> i V r mm (t)l' ram = — 



> 



Q 



xp n {t) 



Thus, using a coupling argument, one can define the processes 
X mm and X min (f) in such a way that X mm (f) < X min (f) for all 
f >0. ■ 

The previous theorem motivates the study of the min policy, 
as it implies for instance that mean transfer times T r under the 
min policy provide upper bounds on the corresponding transfer 
times under the max-min policy. 

III. Ergodicity conditions 

In the following we demonstrate that min and max-min band- 
width sharing policies have a stationary regime under the usual 
conditions, i.e. when the load on each link £ is less than 1: 

Theorem 2: Under the allocation policies £ mm and ^ ) mm , the 
network is 

(0 ergodic if max^ e x Pe < 1 ; 
(if) transient if max^ex pe > 1 . 

This result has already been proven for the max-min policy 
in |6| ; we note that, by Theorem [U ergodicity under the min 
policy implies ergodicity under the max-min policy, thus the 
treatment of the min policy given below provides an alterna- 
tive proof to that of 1 6 1 . However we feel that, since the proof 
below is simpler and uses only elementary Lyapunov functions 
results, it should be easier to adapt to a more complicated sit- 
uation. Transience under condition (ff) is in fact valid for any 
allocation policy which meets the capacity constraints (Q3. 

Proof: Consider the discrete time chain (X(n),n € IN) 
describing the sequence of states visited by the continuous time 
jump process X. Transitions from a given state X = (% r , r € 
satisfy 



JP[Ax r (n) = l\ X(n)=X] - 
W[Ax r {n) = -l \x(n)=x] - 

where Ax r (n) = x r (n + 1) — x r (n) and 



K 

D ' 

1 X r ■ C( 

min — . 

D G r ter Xf ' 



D ^ Y(x r + -mm^-) 

< \H{\- (max A,,- + max maxQ^ 



D'. 



Ergodicity of the continuous time process X will follow from 
that of X and from the fact that the mean sojourn times in each 
state X are bounded from above uniformly in X (or equivalently, 
that the jump rates out of each state X are bounded away from 
zero), a property which is easily verified. 



Sufficient condition. Assume Pm = max^ G £pf < 1, and define 
the Lyapunov function 

/(*)=! E i, 

reH.l<k<x r 

where y r > 1 will be chosen later. The structure of the func- 
tion, which may seem unnatural, has been chosen for the sake 
of computation; it is in fact of the general form EreK. P'Yr r + K, 
for appropriate constants p\ and K. 

In order to express the transition rates in terms of x r and pp, 
remark that 



Xf = Yx r = YX r o r ^— < p^Qmax; 1 
^ "=* A. r O r 



Then, using the notation 

a* <Jcf x r 



r3i X r O r 



X r o r ' 



rent 



we have, Mr G 



JP[Ax r (n) = -1 | X(n) =X]> 



Dp M x 



M 



Thus, 

JE[f(X(n + l))-f(X(n))\x(n)=X] 

= E (y; r+1 P[Ai,-(«) = 1 | X(n) = X] 



reft. 



T fr IP[Ai r («) = -l \X(n)=X]^j 



< E - 

" ,^PmD 



Yr' 



pMjr - TT 
X M 



Let y r = y^57 , where y is such that 

i 

PMjr = PmY 1 ^ < 6 < 1, r e 31, 

for some real number satisfying Pm < 8 < 1. The following 
inequality sums up what we have so far: 

W.[f(X(n+l))-f(X{n))\X(n)=X] 



< Y 

~ f^PuD 



X M 



Let a be a real number such that < a < 1 . The following 
quantities will be evaluated separately 



vi del 

2-2 = 



PmD 



e 



Since Ei is a sum of negative terms, the following bound 
holds, for any ro is such that x* Q = x* M , 



Ei < -5^L_(e-a) < - — (9-a)minA, r < 0. 



PmD 



PmD' 
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Bounding E2 is straightforward: 



^2< E 



6 < -L — ^'^I^IQmaxX,-. 



p M D p M D 



re% 



Now, if C > and e > are chosen to satisfy the inequality 

(e-a)rmnA, r + y (a ~ 1)c |^J9maxA, ; . < -e, 

re%_ re%_ 

we have, MX <E {x* M > C}, 

m[f(X(n + l))-f(X(n))\x(n)=X] 

= Ei +Eo < -e- — < -e— '— < 0. 

p M D p M D' 

Since {x* M < C} is a compact set, Foster's theorem applies 
(see e.g. Q) and the Markov chain is ergodic. 

Necessary condition. Assume now that there exists £0 such that 
pe > 1. Defining 

g (X) = £ o,-x r , 



we immediately have, 

m[g(X(n+l))-g(X(n))\X(n)=X] 

£ 0,-(p[Ai r («) = 1 I X(n) =1] 

-F[Ai,-(«) = -1 I #(«) =X] 



r3(a 



> jy[Ce Pe -Ce }>0. 
Since the jumps are bounded, the chain is transient. ■ 

IV. Mean field analysis of large networks 

It does not appear possible to obtain closed form expressions 
for the stationary distribution of the best-effort network state un- 
der the min policy. We therefore turn to the study of these sta- 
tionary distributions under a limiting regime on network size and 
topology. A similar approach has previously been successfully 
applied to loss networks (see JS), J9), and references therein), 
and to queueing networks in ifTUl , ifTTl . It is inspired by the 
so-called mean field models of statistical physics. 

Mean field analysis in the present context is best illustrated 
by the star-shaped network of Fig. [2] This network has N/2 
branches, each consisting of one inbound and one outbound link 
(thus implicitly N is an even number), and all links have unit 
capacity. Each route r connects the endpoints of two branches 
via the center node. It has an associated arrival rate 2X/N and 
mean message length a. The factor 2/N is introduced to make 
the total load on each link, p = Xa, independent of the number 
of links N. As discussed below, when N goes to infinity the 
number of ongoing transfers on any link becomes independent 
of the number of ongoing transfers on any finite collection of the 
other links (this was termed the "chaos propagation" property 
in ||9])- This allows us to derive fixed point equations for the 
probability distribution of the number of transfers in progress 
on any link. 




connection 



Fig. 2 

A STAR-SHAPED NETWORK 



A. Symmetrical star-shaped networks 

Although amenability to a mean field analysis is a significant 
motivation for considering the star-shaped topology, it should 
be noted that it is also relevant to the study of real networks. 
Any overprovisioned links in a real network are largely trans- 
parent to the throughput of elastic best effort flows. Only bot- 
tleneck resources, typically located in the access network and 
within Web servers, need to be included in the network model. 
The star-shaped network may thus be considered to represent 
any network with a well provisioned backbone where through- 
put is limited by bottlenecks at the source and destination edges. 
For example, inbound links might represent Web server CPU, 
while outbound links correspond to the last hop of an ISP's in- 
terconnection network. This discussion not only motivates the 
consideration of such a topology, but also suggests that letting N 
go to infinity might indeed be realistic if N represents the num- 
ber of Web servers over the Internet. Of course, there would be 
no reason in practice to assume symmetry. This assumption is 
introduced solely for reasons of tractability. 

Although our focus is on the star-shaped topology, the mean 
field approach can be applied to other symmetrical topologies. 
It thus allows one to consider routes with more than two hops. 
The corresponding extended model is described in detail in Sec- 
tion [lV32below, where the corresponding fixed point equations 
are derived. Section lTV-Cl then presents analytical results for the 
star-shaped network. Results of numerical investigation of the 
fixed point equations are reported in SectionlVl 

B. Fixed point equations for large symmetrical networks 

We use the following notation in the sequel. 

• N: total number of links; 

• L: length of a route through the network; 

• R^" 1 : number of routes going through a given link; 

• xf^ : number of active connections on link £ G L, in stationary 
state; 

• Xr : number of active connections on route r, going through 
links r(l),...,r(L); 

• X: arrival rate on a link; 
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Fig. 3 

A HYPERCUBE- SHAPED NETWORK WITH DIMENSION 5. 



• a: mean message length; 

• p = A,o: load of a link. 

We have implicitly assumed here that the number of routes 
going through a link, R [N \ is the same for all links. We shall 
in fact assume further that the network topology is the same, 
as seen from any route. We do not attempt to give a formal 
definition of this symmetry assumption here. The reader is re- 
ferred to [9 1 for a thorough discussion on the minimal symme- 
try assumptions required. Symmetry implies notably that each 
route has the same number of hops and the same traffic pa- 
rameters. The star-shaped network discussed above constitutes 
an example of such a symmetrical network when L = 2 (with 
RW=N/2). 

It is more difficult to come up with meaningful examples of 
symmetrical networks supporting routes with L > 2: in partic- 
ular, the network should not be fully connected, since routes 
longer than one hop would then be pointless. One reasonable 
model is a hypercube (Pig. [3j of large dimension, in which each 
edge contains two one-way links. 

The hypercube is a classical structure with many symmetries. 
It is characterized by its dimension d. Its vertices are represented 
by of-tuples of 0s and Is (e.g., (0, 1, 1,0, 1)) and its edges connect 
two vertices differing in only one coordinate. 

The total number of links in such a network is N = d2 d . The 
number of routes going through any link is 



R w = L 



(d-L)V 



where the only routes considered are the shortest paths between 
two vertices which differ in exactly L coordinates. Note that 
the results below do not depend on the precise topology of the 
network. 

We now derive the fixed point equations. It should be stressed 
that this derivation is heuristic. We clearly mention which steps 
need further justification in the course of the derivation. We do 
believe that the equations are very good approximations, how- 
ever, especially in view of the numerical and simulation results 
presented in the following section. 



Assume now that p < 1 and that the system is in stationary 
state X^ N] = (xf 1 , r G For any k > 0, the proportion of links 
in state k is 



a. 



(N) t* 



By symmetry, it holds that 

P(xf =k) = JEa [ k N \ WeL. 

The chaos propagation assumption^ implies that a[ w) obeys a 
law of large numbers: 



lim ocf 5 = lim JP(Xf> = k) = W(X = k) 



(N) 



(AO 



It appears that the dynamics of the system are driven by OK 
traditionally referred to as the mean field. The following nota 
tion will also be useful: 



a 



(AO 



k>0 



(AO 



a = mx, 



k>0 



In order to derive the equation satisfied by the limit stationary 
distribution OC^, we must first describe the possible transitions 
for ay. The two cases of interest are 

• arrival on a link with k connections: 

a {A0 , „(N) _ 1 (AO 
a k ^ a k iy> U /i+l 

• departure from a link with k > connections: 

1 (AO (AO , 1 



a 



(AO 
a k+l 



1 

+ N 



The transition corresponding to a new connection arrival has 
rate X. The main problem is to compute the departure rate from 
a link i, given that it has xf ] = k ongoing transfers. This can be 
written as 



1 



(A 7 ) 



(4) 



Since the total number of routes is much larger than the size 
of the network, we assume that the probability of having more 
than one connection on a route r is negligibk^], and that the links 
on route r are independent, conditioned on {xf^ = 1 }H. The first 
property allows to rewrite (0]) as 



1 



mm - 

Per X 



(AO 

r(i>) 



4"' = 1 



(5) 



Let j be a given link and let r be one route using link /', i.e.. 



r 3 j. The distribution of Xj N} , conditioned on there being one 
connection on r, is then 



IP[Yj w =Jf |4" } = 1] = 



E[i {x (A0 = 



Ell, (AT) 
{X r - 



We have not proven that this assumption holds. It seems, however, that the 
techniques developed in [ 9] could be applied to prove that this is the case, pro- 
vided the parameter S (A,) goes to infinity with N. 

2 This fact is easy to prove in finite time, but requires more work for the sta- 
tionary regime. 

3 This is the point where the heuristic is not completely exact; it is however 
likely to be true when p tends either to or 1 . 
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By symmetry, it is possible to sum both sides of the above 
fraction over all the routes going through link j, 



TP[xf ) =k'\x^ = l] = 



JEX 
IEoc^' ' 



Departure rate © then becomes, in view of the assumed in- 
dependence property between the X r r{\ given x r = 1, 

Taking the limit AT — > °°, the invariant measure equations fol- 
low. We have: 



A,ao + — oci 



a a 



L-l 







and 



1 / Uk+l Uk 

A, an -a k ) + -( a M —j-^r - o^^— r 



= 



for k > 1 , where 



1 A 1 A 
k k2 



1 



(6) 



(7) 



(8) 



as 



k 2 ,...k L = l "i/ i=2 

Equations © and (|7]l can be rewritten in a more concise form 

Vk > 0. (9) 



V-k+lUk+l — P0C L OCfc, 



The two sets of equations dHJ and (0 together constitute the 
fixed point equations we require. As noted in the introduction 
of this section, these equations do not depend on the topology of 
the network. The expression for Uk can be simplified. 

Let Y2, . . . ,Yi be random variables with distribution 

TP(Yi = y) = ^, 2<i<L, 



Note that, in the case L = 2 (i.e., for the star-shaped network), 
the original equation ([8]l is perhaps simpler than the equivalent 
expression (fTTT >. It yields the following form for the fixed point 
equations: 



V-k+\Uk+l 
Uk 



paa k , 

y>0 



(12) 
(13) 



Remark 1: When considering the star-shaped network as a 
model for Web transfers over the Internet, as suggested in Sec- 
tion IIV-A1 inbound links could be seen as the CPU of Web 
servers and outbound links as the last hop between the ISP's 
backbone and the end customers. It thus makes sense to relax 
the symmetry assumption we had made between inbound and 
outbound links, as the two types of bottlenecks are of a dif- 
ferent nature. We might thus consider a star-shaped network 
with N m inbound links, N M outbound links, inbound (resp. out- 
bound) links having capacity C in (resp. C ml ), see Fig. [4] Assume 
the mean message length a is the same for each two-hop route, 
and the link capacities C°'" are fixed. The arrival rate on each 
route has the form X r = X/N ia , and the load p"'" = A-o/C"" 1 is less 
that 1 . The capacity C in of a "backbone" link is chosen to en- 
sure p in = XoN" M /C n N >n < 1 is fixed when the size of the system 
grows. Then, a&N m andA^ 1 " increase, with N in /N° m small, the in- 
bound links have many active connections and a large capacity, 
while the outbound links remain in "normal" utilization. The 
same approach as above can then be applied, to yield the set of 
fixed point equations 

= P ,n a u, "<, 



where 



u'l = \\ a'° min 

y>0 

U T = £ a v ,min 

y>0 



k,y 



and let Y = max(F 2 , • • -,Y L ). Then © reads 
u k = ka L ~ l IE 



1 ll 

k A Y 



Straightforward calculations yield 



IE 



1 1 

k A Y 



1 ^ 1- 

ytl L £ y- 
1 



(10) 



E 

y=i 

E 



W{Y=y) 

■p(y <y). 



and ai k n (resp. a° k ul ) represents the proportion of inbound (resp. 
outbound) links with k ongoing transfers. 

C. Analytical results for L = 2 

While equation (0 looks superficially like a "birth and death 
process" equation, it is in fact non-linear due to the fact that u\ 
and a both depend on the a v . 

From ([Tol l, one clearly sees that u\ is increasing in k, and 
tends to a L ~' when k — s- 00. Therefore, is increasing as long 
as Uk < pa L ~', and decreasing after that. This means that the oc^ 
form a modal distribution, which maximal value is attained at 



The simplified form for u k is thus, from the basic properties 
of the minimum of independent random variables: 



u k 



k Y — Y ma, 



(ii) 



£0 = max{£ > 0, uk < pa }. 

We now present analytical results on the solution of the fixed 
point equations for L = 2. The proof of these results can be 
found in lfT2l . It relies heavily on functional analysis. 
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Fig. 4 

THE ASYMMETRICAL STAR-SHAPED NETWORK 



Fig. 5 

Distribution of the mean-field probabilities P(X = k) for L = 20 

AND DIFFERENT VALUES OF THE LOAD p. 



Equations (fl~2b — (Tl~3T> have an unique solution for p < 1 and, 
when p — > 1, the following asymptotic expansions hold: 



JEX 



rimp- A 'P(X = A:) 



1 



(1-P) 2 A' 
(1 - p)Bexp 



1 



(l-p)A 



where A and B are non-negative constants. 



Thus, under the min policy, any link in the star-shaped net- 
work has a mean queue length which is one order of magni- 
tude larger than for a single server queue with the same load 
(p/(l — p)). Its tail distribution is still geometrical with factor 
P- 

It is possible to give an expression for the constant A: if c and 
v are solutions of the following system of differential equations, 



zc'(z) + c(z)v(z) = 0, 
zv"{z) + v'(z)=c(z), 
[v(0) = 0, v'(0) = l, 



c(0) = l, 



then A can be written as follows: 



c{z)dz= lim zv'(z) « 1.30. 



Since this system is numerically highly unstable, it has proven 
difficult (with the "Livermore stiff ODE" solver from MAPLE) 
to derive a better estimate for A. 

It is interesting to note that the function w(y) = v(e y ) + 1 sat- 
isfies the so-called Blasius [13] equation 

w"'{y) + w{y)w"{y) =0, 

which describes a laminar boundary layer along a flat plate (see, 

e.g., mi). 

V. Numerical analysis and simulations 



While the analytical results of Section HV-CI give some good 
estimates, they are only valid in the heavy traffic regime p — > 1 . 



In addition, similar results for L > 2 are not available. We thus 
resort to numerical resolution of the equations to gain a better 
understanding of the performance of transfers across large sym- 
metrical networks. The very form of the equations suggests the 
use of a fixed-point method for this numerical resolution: start- 
ing from a priori values (a[ 0) , k > 0), the algorithm computes the 
corresponding m[ 0) from ®, and then new values a^ 1 from (|9). 
Provided special care is taken to avoid instabilities, the iteration 
of this process converges rapidly (less than 100 steps). Sample 
results are shown in Fig. [5] corresponding to a large symmetrical 
network with routes of length L = 20. 

As clearly seen in the figure, the distributions are very dif- 
ferent from what would be obtained for routes of length L = 1 . 
In this case, the system consists of a collection of independent 
M/M/l queues and the associated distribution {o^} is geometric. 
For L > 1 , the {oc^} distributions are markedly modal and the po- 
sitions of the peak values are roughly proportional to (1 — p)~ 2 , 
a fact has only been proven in Section ITV-CI for the case L = 2. 
Moreover, since the shape of the distribution is rather narrow, 
this position roughly coincides with the mean number of active 
connection (as can be seen from the raw data). 

The impact of route length is illustrated in Fig. [6] It seems 
that the mean number of active connections (which is again ap- 
proximately the peak value of the distribution) is roughly pro- 
portional to logL. Note that a logarithmic growth rate is very 
slow suggesting that, beyond 2 or 3, the number of bottlenecks 
does not have a significant impact on mean transfer times. They 
depend much more on the load p. 

The results presented so far only concern the solution of the 
fixed point equations. As mentioned earlier, there are gaps in 
the derivation of these equations. To assess their quality and 
to investigate the accuracy of the asymptotic approximation for 
finite size networks, we ran a number of simulations of the star- 
shaped network. Fig.|7]displays the corresponding results when 
the load on each link is set to p = 0.9 for a varying number 
of links. The agreement between the simulation results and the 
fixed point equation results is excellent for N = 100 links and 
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Fig. 6 

Distribution of the mean-field probabilities F(Z = k) for 

DIFFERENT VALUES OF THE ROUTE LENGTH L AND LOAD p = 0.9. 
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Fig. 7 

Cumulative distribution F(X < k) for different values of N 

(SIMULATION) AND INFINITE SIZE (FIXED POINT) FOR A STAR-SHAPED 
NETWORK WITH LOAD p = 0.9. 



improves as N increases. 

VI. Conclusions 

In this paper, we have considered a class of Markov pro- 
cesses called best-effort networks which constitutes a natural 
probabilistic model for evaluating the performance of document 
transfers over data networks such as the Internet. Unlike almost 
all previous work, this model accounts for the random nature 
of traffic: document transfers begin at the epochs of a certain 
arrival process and the size of each document is drawn from 
a given probability distribution. In the interests of tractabil- 
ity we assumed Poisson arrivals and exponentially distributed 
sizes. We introduced the "min" bandwidth sharing policy as a 
conservative approximation to the more classical max-min pol- 
icy. Necessary and sufficient ergodicity conditions for best ef- 
fort networks under the min and max-min policies have been 
established. 



In order to pursue the analysis of the stationary distributions 
of the number of transfers in progress, we have resorted to large 
network asymptotics applying the mean field approach of sta- 
tistical physics. This enabled us to derive fixed point equations 
for the probability distribution of the number of ongoing trans- 
fers on a given network link. The validity of these equations has 
been established by comparing their solution with the results of 
simulations. 

Analytical and numerical results show how the mean transfer 
time depends on the number of bottleneck links and their load. 
The steady state distribution in networks where routes have sev- 
eral bottlenecks (L > 1) has a marked modal behavior. This is 
significantly different to the geometric distribution which holds 
when routes have a single bottleneck (L = 1). Performance is 
also much more sensitive to link load p for multiple bottleneck 
routes: as p — > 1, mean transfer time increases like 1/(1 — p) 2 
in the case L = 2, whereas the dependence is in 1/(1 — p) when 
L = 1 . Finally, the impact of the number of hops per route L 
appears small (given that L > 1) compared to that of parame- 
ter p. This suggests that the star-shaped network is perhaps a 
sufficiently complex model, and that the study of symmetrical 
networks with L > 2 is less relevant. 

The work presented here can be pursued in several directions. 
On the theoretical side, the analytical results presented in Sec- 
tion IIV-CI constitute a first step to understanding the solution 
of the fixed point equations which could be taken further. An- 
other challenging theoretical question is to improve the fixed 
point equations in a rigorous way. On a more practical side, the 
fixed point equations might be simplified so as to find simple 
approximate formulas for mean transfer times as a function of 
key parameters (such as p'", p oul in the case of the asymmetrical 
star-shaped network described in RemarkQ}. Such approximate 
formulas could then lead to engineering rules for capacity plan- 
ning. 

We view the present study as a preliminary investigation into 
the performance of best effort networks with multiple bottleneck 
links. A significant result of this investigation is the discovery 
that the extension of the processor sharing model valid for a 
single bottleneck proves to be very hard. There appears to be 
no simple parallel to the familiar fixed point techniques used in 
loss networks. The problem is, however, of considerable prac- 
tical importance for providers seeking to engineer their network 
to ensure adequate throughput for document transfers. We hope 
therefore that this paper will incite further work and the devel- 
opment of alternative heuristic approaches. 
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